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Abstract. Digitized images of full disk CaK spectroheliograms from two solar observa- 
tories were used to study cycle variation of ephemeral regions (ERs) over ten solar cycles 
14-23. We calculate monthly averaged unsigned latitude of ERs and compare it with an- 
nual sunspot number. We find that average latitude of ERs can be used as a predictor for 
strength of solar cycle. For a short-term prediction (dT ~ 1-2 years), maximum latitude of 
ephemeral regions (in current cycle) defines the amplitude of that cycle (higher is the lati- 
tude of ERs, larger are the amplitudes of sunspot cycle). For a long-term prediction (dT ~ 
1.5 solar cycles), latitude of ERs at declining phase of n'' 1 cycle determines the amplitude of 
(n+2) f ' ! sunspot cycle (lower is the latitude of ERs, stronger is the cycle). Using this latter 
dependency, we forecast the amplitude of sunspot cycle 24 at W=92±13 (in units of annual 
sunspot number). 
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1. Introduction 

Historic data sets of full disk Ca II K 
spectroheliograms observed from three obser- 
vatories: Kodaikanal (KKL), Mount Wilson 
(MWO), and the National Solar Observatory 
at Sacramento Peak (NSO/SP) span about ten 
past solar cycles. Recently, these data were 
digitized and calibrated. In this study we use 
KKL and NSO/SP datasets to explore latitu- 
dinal distribution of ephemeral active regions 
as potential precursor of amplitude of sunspot 
cycles. Images were acquired with exit slit 
of spectroheliograph of about 0.5 A in width 
centered at A = 3933. 67A. Spatial resolution 
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is approximately 1.2 arc seconds per pixel. 
Fu rther details on these data sets can be found 
in iTlatov et all d2009l) . Flocculi and plages 
were ident ified using intensit y threshold as de- 
scribed in lTlatov et al.l d2009h . Isolated clusters 
of bright pixels (not connected with other clus- 
ters) were classified as independent flux ele- 
ments, and the total number of elements and 
area of each element were computed for ev- 
ery image in our data set. Despite differences 
in observations, floccular areas show good cor- 
relation between ins truments and with s unspot 
number (Figure Q}. Tlat ov et"ai] ([2009) estab- 
lished a correlation between total area of CaK 
plages (Api age ) and area of sunspots (A slmspot ), 
Apiage = (8.5 ± 0.3) + (15 + 0.25)A SB „ spor , with 
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Fig. 1. Monthly averaged area of Ca II K plages de- 
rived from data from three observatories (three up- 
per panels) and area of sunspots (lower panel). 

Pearson linear correlation coefficient R=0.88. 
KKL data set is the longest of three, while 
NSO/SP set is the most recent. Given good cor- 
relation between KKL, MWO, NSO/SP Ca II 
K data we chose to use the Kodaikanal and 
NSO/SP data for the following analysis. 

2. Average Latitude of Ca II K Flocculi 

We have selected a subset of mid-size floc- 
cular elements with area between 100 and 
300 millionths of solar hemisphere (MSH). 
lHarvev & Martin I d 19731) had classified ele- 
ments with this area observed in Ca II K line as 
ephemeral regions. Majority of bright elements 
identified in Ca II K images correspond to the 
chromospheric network elements; only a small 
fraction (~ 10%) matches ephemeral regions 
and plages. Furthermore, while total number of 
bright elements does not vary with sunspot cy- 
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Fig. 2. Number of Ca II K ephemeral regions in 
NSO/SP data as function of time and latitude. Data 
corresponds to 3-monfhs averages computed over 5° 
latitudinal intervals. 



cle, number of ephemeral regions does exhibit 
sunspot cycle-like variations. Figure [2] shows 
number of ERs found in NSO/SP data as func- 
tion of time and latitude. Typically, cycles of 
ephemeral regions begin at latitudes of ~ 60 
degrees about one-two years before sunspot 
maximum. ERs cycles appear to be shifted rel- 
ative to sunspot cycles: number of ERs is at its 
lowest at the rising phase of sunspot cycle, and 
it reaches maximum in the tail of sunspot cycle 
when sunspot activity concentrates at the equa- 
torial region. 

Next, we have calculated average latitude 

N 

of ephemeral regions 6 — jj ^ where N is 

total number of ERs observed in a given month 
and f?, is the latitude of individual ER. 

Average latitude of ephemeral regions 
shows clear variation with sunspot cycle 
(Figure [3}. ERs latitudes reach maxima on the 
rising phase of sunspot cycle, approximately 
one - two years before the number of sunspots 
reaches its maximum. Examining Figures [3] 
and|4]one may note a correlation between max- 
imum latitude f?i of ephemeral regions and am- 
plitude of sunspot cycle: higher is the average 
latitude of ERs, stronger is the sunspot cycle. 
Similar tendency is also present in the butterfly 
diagram of sunspots (not shown): in stronger 
cycles first sunspots emerge at higher latitudes, 
while for weaker cycles, latitude at which early 
sunspots emerge is lower. Figure [5] shows cor- 
relation between maximum latitude of ERs {Q\ ) 
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Fig. 3. Three months (upper panel, thin curve) and 
12-monfhs (thick curve) averaged latitude of Ca II K 
ephemeral regions derived from NSO/SP data and 
annual sunspot number (lower panel). Thin verti- 
cal lines demonstrate that ephemeral regions reach 
maximum latitude at rising phase of sunspot cycle, 
a few years prior to sunspot maximum. 



and sunspot number (W); Pearson linear corre- 
lation coefficient Rgi = 0.83. 

Even more interesting, minima of average 
latitude of ERs (62) that occur on declining 
phase of n' h cycle correlate quite well with am- 
plitude of (n+2)'' 1 sunspot cycle. Dashed lines 
with arrowheads shown in Figure [4] demon- 
strate this relation. 

Figure [6] shows correlation between mini- 
mum latitude of ERs (62) and sunspot number 
(W) with Re2=0.92. This latter correlation en- 
ables forecast of amplitude of sunspot cycle for 
about one and a half solar cycles prior to its oc- 
currence. 

3. Discussion 

Forecasting the strength of future solar cy- 
cles was attempted by number of r esearchers 
(e.g., lor] 119681: LTavaraiahl 120071) . In our 
present study, we demonstrate that average 
latitude of ephemeral active regions can be 
used for a short-term (dTi ~ 1-2 years) and 
a long-term (dT2 ~ 1.5 solar cycles, or 14-17 
years) forecast of amplitude of sunspot cycle. 
Appearance of emerging regions at high lat- 




1980 2000 
Years 



Fig. 4. Average latitude (9, in degrees) of ephemeral 
regions (upper panel) and sunspot number (W, lower 
panel). Open squares and circles mark local min- 
imum (ft,) and maximum 0ij) latitudes of ERs. 
Dashed lines with arrowheads draw a correspon- 
dence between minimum latitude of ERs and max- 
imum sunspot number. Solid lines with arrowheads 
show correspondence between maximum latitude of 
ERs and maximum sunspot number. To see the pat- 
tern, compare order of cycles in their amplitude with 
order of minima and maxima in average latitude of 
ERs. 
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Fig. 5. Amplitude of n''' sunspot cycle (annual 
sunspot number) as function of maximum latitude 
(in degrees) of ERs in rising phase of same cycle. 
Solid line shows a linear approximation to the data, 
W" = (-445.782 ± 146.000) + (24.884 ± 6.289) x ff[ . 
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Fig. 6. Amplitude of n' A sunspot cycle (annual 
sunspot number) as function of minimum latitude 
(in degrees) of ERs in (n-2)''' cycle. Solid line shows 
linear approximation to the data W" = (541.467 ± 
69.939) - (26.761 + 4.622) x e:p 2 . 

itudes heralding emergence of sunspots is in 
agre ement with the idea of extended solar cy- 
cle (IWilson et al.l ll988) with magnetic activ- 
ity of a current cycle starting at high latitudes 
a few years prior to sunspot emergence. In its 
turn, dT2 is very close to a time interval to 
transport surface magnetic field into the dy- 
namo region a s required by the transport dy- 
namo models fchoudhur i et al. Ifl995b iTlatov I 
1996). Latitude 62 is lower, when the old cy- 
cle activity penetrates to lower latitudes in ab- 
sence of high-latitude activity of a new cycle. 
This implies that 62 is lower, when overlap be- 
tween old and new cycles is shorter. Since the 
orientation of magnetic fields in two consec- 
utive cycles (cycles n and n+1) is opposite to 
each other, longer overlap between cycles will 
lead to a weakening of a seed field for (n+2) 
cycle. In addition, a penetration of magnetic 
field to lower latitudes may allow for a more 
efficient cancellation of magnetic field of oppo- 
site polarity across solar equator, which in the 
framework of Babcock-Leighton model may 
lead to a more effective transformation of mag- 
netic fields of active regions into poloidal field 
of future cycles. 

Finally, using found relationship between 
minimum latitude of ephemeral regions and 
amplitude of solar cycle (Figure |6j, we pre- 
dict a relatively modest cycle 24 with annual 
averaged amplitude of about W=92+13. This 



our prediction is in agreement with recent fore- 
cast of stre n gth o f Cycle 24 (92.8+19.6) by 
iBhatt et a!71 d2009l) based on geomagnetic aa 
index. 
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